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1. INTRODUCTION 

Recently, the all-electric ship concept has become particularly important. Full electric propulsion 
systems for commercial and military shipping sector have been intensively developed. Electric propulsion 
offers potential advantages regarding flexibility in ship maneuvering, maintenance cost and vibration level 
[1], [2]. In the naval industry, three-phase electric machines are mainly the most used. Recently, the use of 
multiphase machines has grown considerably in several areas including naval propulsion, railway electric 
traction, electric vehicles, avionics, and high-power industrial applications [3]-[9]. Multiphase machines 
provide several important benefits in terms of torque quality and machine reliability [1]. The most interesting 
multiphase machines which attract a lot of interest is the dual three phase (DTP) ones. DTP machines has two 
sets of three-phase stator windings spatially shifted by 30 electrical degrees with two isolated neutrals [1]. 

Among all DTP machines the permanent magnet synchronous motor (DTP-PMSM) is the most used 
one [10]. A field oriented control of DTP-PMSM based on the vector space decomposition (VSD) was 
presented in [10]. The conventional VSD control for DTP induction motor was given in [11]. As specified by 
the VSD strategy, the machine model is transformed into three decoupled subspaces, identified as, (a, p) 
torque-component, (Z1, Z2) harmonic-component and (01, 02) zero-sequence, respectively [12]-[18]. 
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Direct torque control (DTC) is a technique used for high performance control of three phase electric 
drive systems [19]. It is characterized by simple structure which is independent of machine parameters and 
allowing a fast torque response. When applying the DTC for the DTP-PMSM, important harmonic stator 
currents are usually observed. These currents cause losses in the stator and thus deteriorate the machine’s 
efficiency. According to the VSD technique, the basic DTC does not allow the control of the harmonics that 
appear in the subspace (Z1, Z2). 

With the view to minimize harmonics and hence improve the machine’s efficiency, a modified 
Switching Table based DTC for five phase synchronous machine is proposed in reference [20]. In this 
approach the selection of the appropriate voltage is done in two steps. The same technique applied to DTP- 
PMSM is presented in reference [21]. The disadvantage of this method is the necessity to locate the stator 
flux position not only in the (æ, f) subspace but also in the (z1, z2) subspace. Our contribution is to propose a 
modified switching table constituted by synthetic vectors allowing the torque control as well as the reduction 
of the currents in (Z1, Z2) subspace. 


2. DTP-PMSM DYNAMIC MODEL 

Figure 1 and Figure 2 illustrate the DTP-PMSM and VSI-fed drive. The two three-phase windings 
are identical with two independent neutral points [12]. The DTP-PMSM is complex and high order system. 
The mathematical model of the motor is simplified by assuming that the motor windings are distributed 
sinusoidally, saturation and magnetic losses are neglected [22]. To obtain a practical model proper for 
control, the voltage space decomposition VSD is used [11]. The VSD is based on the transformation matrix 
expressed by (1). 





Figure 1. DTP-PMSM ship and six-phase voltage Figure 2. Windings of DTP-PMSM 
source inerter 


1 V3 V3 
Po a SS SC 
V3 V3 1 1 
oe 2 ea 
1 4 1 V3 v3 
U “3 “4 = = » (1) 
y B 1 4 i 
2 2 2 2 
1 1 1 0 0 0 
0 0 0 1 1 1 


Using this transformation, the system is projected into three orthogonal subspaces (a, P), (Z1, Z2) and 
(01, 02). The voltage equations, are expressed as [12]: 


[Vag] = [Rs] [iag] + < [Wap l=[Rs] [iag] + < [[Lag] [ieg] + Wpm- i A (2) 
[Vz] = [Rs] liz. | + < lw, | — [Rs] liz, | + [Lz] < [iz] (3) 
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d dy. 
aa = [Rs] Eas + dt [Poz] = [Rs] lio. a| + [Lo] dt [iona] (4) 
Where: 
Se Sp e 5596 ee sin 20 
[tap] = seu sin 20 Ae — Cato cos 26 


La, Lg are d and q axis inductances. 
L; and L, are the transformed stator self-leakage inductance. 
W py 1s the PM flux linkage. 
0 is the rotor position. 
Fundamental and harmonics are distributed as; i) the fundamental component is projected in (a, B), 
ii) harmonics with order 12K*" + 1(K = 1, 2,3,...) are projected in (a, B), iii) harmonics with order 6K" + 
1 (K = 1,3,5,...) are projected in (zı, z2), and iv) the harmonics with order 3K" (K = 1,3,5,...) are 
projected in zero sequence subspace (01, 02). 

Current components in the (a, p) subspace produce air gap flux and create the torque. Currents in 
(Z1, Z2) and (01, 02) are all harmonics, they don’t contribute to the electromechanical conversion and produce 
stator losses [23], [24]. To change to the (d, q) rotating frame the following transformation matrix is applied: 


_ [cos sind 

r = L-sin@ Be ©) 
Machine model in the (d, q) plane: 

Va _ YR; 0 lg a [Wa d0 Pa 

vsl=[0 rollig talba La. © 

Wa oe La 0 Íg Wpm 

wl=lo if) tele" (7) 
The electromagnetic torque is: 

Te = p(igWa = iapa) (8) 


p : pole pairs number. 


3. RESEARCH METHOD 
3.1. Inverter voltage vectors 

The inverter shown in Figure 1 can provide 64 different voltage vectors. Each vector is presented by 
the subscript decimal number, corresponding to binary numbers SaiSpiSciSa2Sp2Sc2. Where (S = Sai, Spi, Sci, 
Sa2, So2, Sc2) are the device switch states. The phase voltages are expressed according to switch states by (9). 


Va 2 -1 -1 o0 0 0 qa 
V1 —1 2 -1 0 0 0 ||% 
Vaj _Eļ-1 -1 2 0 0 0 ||Sa (9) 
Vol 3| 0 0 0 2 -1 —11/S,5 
Vio 0 0 0 -1 2 -1]|S,5 
ae 0 0 0 -i -1 2415, 


While (S = Sai, Spi, Sci, Sa2, Sb2, Sc2) are the states of switches, E is the DC voltage. The 64 voltage vectors 
projected in each subspace are given by (10). 


[Ve Vg Ver Veo Vor Voz] =[T][Var Vor Ver Vaz Vo2 Veo]? (10) 


Thus, in (a, J) and (z1, z2) subspaces, there are sixty non-zero voltage vectors and four zero voltage vectors 
(Vo, V7, Vz6, V62) [25] as shown in Figure 3 and Figure 4. 
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Figure 3. Voltage vectors diagram in (a, f) Figure 4. Voltage vectors diagram in (Z4, Z2) 


According to the Figure 3, the voltage vectors are distributed in the (a, /) plane forming four dodecagons 
(D,, Dz, D3, D4) [26]. Their amplitudes are given by (11). 





(2-v3) 
1 
Up2 = 3 
2 (11) 
Up3 = V3 
(2+V3) 
Ups = 


We can notice from Figure 3 and Figure 4 that vectors of the outer dodecagon D; in (a, p) subspace 
map in the inner dodecagon D; of the (z1, z2) subspace, the inner dodecagon D; of (a, p) subspace forms the 
outer dodecagon Ds, of (z1, Z2) subspace while the middle dodecagon space vector keeps the the same region. 


3.2. Conventional ST-DTC 

In the direct torque control scheme, torque and flux of DTP-PMSM< are estimated. Two hysteresis 
comparators for torque and flux and a switching table are used to generate the inverter vectors as shown in 
Figure 5. The voltage vectors with the maximum amplitude divide the subspace (a, p) into 12 sectors as 
shown in Figure 6. This choice allows maximum use of the DC power supply and guarantees the application 
of low amplitude voltage vectors in the (Zi, Z2) subspace. Each sector is delimited by two maximum vectors, 
as shown in Figure 6. When the stator flux is in the sector k, the application of the vector Vķ+2 leads to the 
increase of the flux and the torque whereas the application of Vg-4 causes their reduction. The application of 
the vector V,43 produces the reduction of the flux and the increase of the torque, the opposite effect is 
produced by the vector V,_3. In summary, the choice of the voltage vector is carried out according to the 
Table 1. 
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Figure 5. Conventional DTC diagram of DTP-PMSM drive 
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Figure 6. 12 sectors diagram and selection of voltage vectors 


Table 1. Conventional DTC switching table 


K Sector 
Hy 1 i 
Hr, l 0 -1 l 0 -1 
Applied 
ne Vk+2 Vero Vk-3 Vk+3 Wera Vk-4 


The generated torque and flux control signals Hy, and Hr, are defined as: 


lif T — T, = Er 

Hr, =} Oif T-T, =0 
aiT- 
raion ee 

—1lif Ys — W, > Ey 


4. THE MODIFIED ST-DTC APPROACH 


The modified switching table DTC approach allows the decreasing of the current harmonics by 
reducing the currents in subspace (zı, z2). As shown in Figure 7 the sixty non-zero voltage vectors of the 
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subspace (a, p) can be divided into twelve groups of vectors {G,, G> ..., G42}. Each group is composed of 
three vectors having the same direction, for the example V3, , V;3 and V46 constitute the G; group. Given this 
remark, the vectors of the same group will have the same effects on the flux and the torque. 

The vectors belonging to the same group in (a, P) subspace change the direction and the module in 
(Z1, Z2) subspace. Indeed, for a given group the largest vector of the Dodecagon D4 changes direction in 
subspace (Zi, Z2) and its module becomes the smallest. The second vector of the group belonging to the 
Dodecagon D3 keeps the same module but changes direction and becomes opposite to the preceding vector. 
As illustrated in Figure 8, the vectors V3, Vig and V3, of D3, Dı and D3 respectively and which constitute the 
group G; have the same direction. However, in the (z1, Z2), V53 has an opposite direction to V3_ and V46 as 
shown in Figure 8. 





Figure 7. Groups of vectors {G,, G3 ..., G12 }in (a, P) Figure 8. Vectors V36 , Vig and Ve3 in (Z4, Z2) 


Consequently, the use of vector of D3 and vector of D4 together allows controlling current in (Z4, Z2) 
and then reducing current harmonics. For the example if V36 is chosen according to the switching table, the 
vector V- must be applied with V3, to cancel its effect in (Z4, Z2) subspace. The (12) relation gives the time 
allotted for each vector: 


lone F T2(Vs53) 2125 = 0 (12) 
T, + T, =T; 


Where T, is the sampling time. The calculation of the times gives: 
Ti = (V3-1)-T, 
(13) 
Ta = (2-V3)-T, 


Thus, the proposed approach consists of replacing the 12 vectors of D4 used in the conventional DTC by a 
combination of vectors of D3 and Dz as depicted in Table 2. The module of synthetic vector is given by (14). 


T1V36+T2Vs53 
TS 


IVI] = = (V6 —- V2) E (14) 














The synthetic vector has a module slightly lower than that generated by the classical DTC. 





— lL X 100% = 92.82% (15) 
IlVsell 
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Table 2. Synthetic vectors used in modified ST-DTC 


Classical DTC Modified DTC 
V36 T1V36 + T2V53 
Vs2 T,Vs2 + T2V3g 
Vz4 Ti V54 + T2V20 
V22 Ti V22 + T2V50 
Vig TyVig F T2 V30 
V26 Ti V26 + T2Vi9 
Vo7 Ti V27 + T2Vio 
Via Ti Vii + T2V25 
Vo TiVo + T2 V33 
Vya Ti Vz F T3 Vis 
Vas TVs + T2V33 
V37 Ti V37 + Tz V4 


o 


1321 


5. RESULTS AND DISCUSSION 

The simulations via MATLAB/Simulink environment were performed using information found in 
Table 3 [21]. The two studied strategies conventional and proposed ST-DTC for DTP-PMSM drives, have 
been tested and results are presented in Figures 9 to 15. As shown in Figure 9, using the classic DTC strategy 
results in a non-sinusoidal phase current, indeed, the current contains a large quantity of the 5°” and 7°” 
harmonics which are dominant (THD = 51.6%) as shown in Figure 11. These current components do not 
contribute to the electromechanical conversion and only generate losses. 

Figure 10 shows the phase current corresponding to the proposed ST-DTC strategy. It can be 
observed that the shape of the current is nearly sinusoidal, harmonic analysis presented in Figure 12 shows 
that the harmonics are significantly reduced (THD about 11%), this is due to the fact that modified DTC 
method allows the control of current components in the (Z1, Z2) subspace, and this consequently allows the 
reduction of harmonic currents. The Figure 13 and Figure 14 represents the currents in (a, p) subspace, we 
notice that for the two approaches we have the same amplitude, low ripple and regular trajectory. In Figure 
15 (a) and Figure 15 (b) the torque presents the same appearance for both strategies. The torque quality is not 
affected by the presence of aforementioned harmonics. We also notice that the velocity reaches its setpoint 
with good dynamic and zero static error. 


Table 3. Simulation parameters 


Designation Value 
DC voltage E 50V 
Resistance R, 1.09Q 
D axis inductance Lg 2.142mH 
Q axis inductance L, 2.142mH 
Permanent magnet flux Wpy 0.0734Wb 
Inertia moment J 89.10°kgm? 
Viscous friction coefficient f 0.01 Nms/rad 
Pole pairs number p 5 
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sN 5 
x 
wa’ 0 
wb 
‘oo 
-10 
-15 
-20 
-25 
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Figure 9. Stator phase current conventional DTC 
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Figure 10. Stator phase current modified DTC 
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Figure 11. Stator phase current harmonic analysis conventional DTC 
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Figure 12. Stator phase current harmonic analysis modified DTC 
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Figure 13. Currents in (a, f) and (z1, z2) subspaces 


Figure 14. Currents in (a, J) and (z1, z2) subspaces 
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Figure 15. Speed and torque dynamic for (a) conventional DTC; (b) modified DTC 


6. CONCLUSION 


DTC is one of the strategies used for high performance electrical drive systems. Applied to DTP- 
PMSM the conventional DTC leads to significant harmonic currents. To deal with this problem, a DTC with 
a modified switching table has been proposed to control the DTP-PMSM. The purpose of this method is to 
minimize the harmonics of the stator current. It is based on the use of 12 voltage vectors synthesized in a 
two-step process. This technique allows the elaborating of the most suitable inverter voltage vector which 
permits not only the control of variables in (æ, p) subspace but also the reduction of currents in the (Z1, Z2) 
subspace. Simulations have shown the effectiveness of the modified strategy to minimize harmonic current 
and increase system efficiency, while preserving the benefits and the merits of the classical technique. 
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